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PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  provides  a 
method  to  calculate  damage  progression  on  a  rubble -mound  breakwater,  revetment,  or  jetty  trunk 
annor  layer.  The  methods  apply  to  uniform-sized  armor  stone  (0.75  W50  <  Wso<  1 .25  (Tso,  W50 
=  median  weight  of  armor  stone)  as  well  as  riprap  (OA25W5o  <  W50  <  4W5o)  exposed  to  depth- 
limited  wave  conditions. 

The  equations  discussed  herein  are  primarily  intended  to  be  used  as  part  of  a  life-cycle  analysis, 
to  predict  the  damage  for  a  series  of  storms  throughout  the  lifetime  of  the  structure.  This  life- 
cycle  analysis  including  damage  prediction  allows  engineers  to  balance  initial  cost  with  expected 
maintenance  costs  in  order  to  reduce  the  overall  cost  of  the  structure.  The  equations  are  intended 
to  provide  a  tool  for  accurate  damage  estimates  in  order  to  reduce  the  possibility  of  unexpected 
maintenance  costs. 

INTRODUCTION:  Rubble-mound  breakwater,  revetment,  and  jetty  projects  require  accurate 
damage  prediction  as  part  of  life-cycle  analyses.  However,  few  studies  have  been  conducted  to 
determine  damage  progression  on  stone-armor  layers  for  variable  wave  conditions  over  the  life 
of  a  structure.  Previous  annor  stability  lab  studies  were  intended  to  determine  damage  for  the 
peak  of  a  design  stonn.  As  such,  most  previous  laboratory  studies  were  begun  with  an 
undamaged  structure  and  damage  measured  for  a  single  design  wave  condition  (Hudson  1959; 
Van  der  Meer  1988).  The  empirical  equations  derived  from  these  studies  were  valid  for  deter¬ 
mining  initial  damage  but  not  for  damage  progression  through  several  storm  events.  Damage 
actually  occurs  as  a  result  of  a  sequence  of  storms  of  varying  severity  and  with  varying  water 
levels.  This  CHETN  provides  equations  that  allow  the  prediction  of  rubble-mound  deterioration 
with  time.  These  relations  are  supplemented  by  predictive  equations  for  the  uncertainty  or 
variability  of  damage  for  more  accurate  estimation  of  reliability  or,  conversely,  probability  of 
failure. 

Within  this  CHETN,  damage  is  defined  in  terms  of  the  average  normalized  cross-sectional 
eroded  area  of  armor  on  the  slope.  Damage  is  defined  up  to  the  point  that  the  underlayer  is 
exposed  through  a  hole  the  size  of  a  nominal  armor  stone  diameter  Dn5 0  =  (M5 o/pa)  ,  where  M50 
is  the  median  mass  of  annor  stone  and  pa  is  the  armor  stone  density.  The  condition  where  the 
underlayer  is  exposed  defines  failure  of  the  armor  layer  because  rapid  destruction  of  the  structure 
often  occurs  after  this  point.  The  damaged  profile  is  described  in  terms  of  the  engineering 
parameters  maximum  eroded  depth,  minimum  remaining  cover  depth,  and  maximum  cross-shore 
length  of  the  eroded  region.  Relations  for  these  profile  descriptors  are  given  in  terms  of  the  mean 
damage.  Further,  relations  describing  the  alongshore  variability  of  damage  and  the  profile 
descriptors  are  provided  to  support  reliability  or  uncertainty  analyses.  These  relations  apply  for 
single  storms  and  for  storm  sequences  given  depth-limited  normally  incident  waves.  The 
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relations  and  supporting  studies  are  described  in  a  series  of  publications  on  damage  (Melby 
1999;Melby  andKobayashi  1998a,  1998b,  1999). 

DAMAGE  DESCRIPTION:  The  Shore  protection  manual  (1984)  provides  damage  as  a 
function  of  the  marginal  wave  height  exceeding  the  zero  damage  wave  height.  The  Shore 
protection  manual  damage  D%  was  defined  as  the  nonnalized  eroded  volume  in  the  active 
region,  extending  from  the  middle  of  the  breakwater  crest  down  to  one  wave  height  below  the 
still-water  level.  This  design  information  is  based  on  laboratory  tests  limited  to  monochromatic 
nonbreaking  waves  impinging  on  long  structure  slopes.  The  background  reports  provide  little 
insight  and  no  data. 

Broderick  and  Ahrens  (1982)  provided  a  definition  of  damage  that  was  not  a  function  of  cross- 
sectional  geometry.  They  defined  damage  to  an  armor  layer  by  the  normalized  eroded  cross- 
sectional  area 


5  = 


(1) 


where 


Ae  =  measured  eroded  cross-sectional  area  (Figures  1  and  2) 

Dn 50  =  nominal  armor  stone  diameter 

This  damage  formulation  was  popularized  by  Van  der  Meer  (1988).  Melby  and  Kobayashi 
(1998a)  utilized  S  as  a  general  damage  description  and  further  defined  the  eroded  profile  using 
several  engineering  parameters:  the  maximum  eroded  depth  E  =  deIDn 50,  the  minimum  remaining 
cover  depth  C  =  dJDn5 0,  and  the  maximum  cross-shore  length  of  the  eroded  region  L  =  le/Dn50. 
These  damaged  section  descriptors  are  shown  in  Figure  1.  A  typical  structure  cross  section  is 
shown  in  Figure  2.  The  maximum  eroded  depth  and  minimum  remaining  cover  depth  may  not 
occur  at  the  same  location  along  a  structure  because  of  the  variability  of  the  armor  layer 
thickness. 


PREDICTIVE  RELATIONS:  Melby  and  Kobayashi  (1998a)  conducted  a  series  of  experiments 
measuring  the  erosion  of  a  stone-armor  layer  for  varying  wave  and  water-level  conditions.  The 
structure  profile  was  measured  repeatedly  throughout  the  test  at  up  to  32  sections  alongshore. 
The  32  profiles  were  used  to  obtain  mean  damage  and  mean  damaged  profile  as  well  as  the 
alongshore  variability  of  damage  and  the  profile. 

The  empirical  equation  proposed  by  Melby  and  Kobayashi  (1998a)  for  predicting  the  temporal 
progression  of  mean  eroded  area  as  a  function  of  time  domain  wave  statistics  is: 


S(t)  =  S(tn)  +  0.025 


(Ns)l 

(T  )0'25 

v  m  /  n 


(t0-25-O  for  tn<t<t, 


.0.25 


n+ 1 
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Figure  2.  Rubble-mound  structure  cross  section 


where  S(t)  and  S(tn)  are  predicted  and  known  mean  eroded  areas  at  times  t  and  tn,  respectively, 
with  t  >  t„.  Ns  =  Hs  /  (A D,75o)  is  the  stability  number  based  on  the  average  of  the  highest  one -third 
wave  heights  from  a  zero-upcrossing  analysis,  A  =  Sr  -1  where  Sr  is  the  armor  stone  specific 
gravity,  and  Tm  is  the  mean  period.  The  wave  parameters  are  defined  5 Hs  seaward  of  the 
structure  toe,  which  is  the  travel  distance  of  large  breaking  waves.  Equation  2  provides  a  means 
to  compute  damage  over  a  sequence  of  N  events,  each  of  constant  wave  conditions,  where  each 
event  is  defined  over  a  time  period  from  t„  to  tn+ 1,  1  <n<  N. 
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A  similar  equation  relating  mean  damage  to  spectral  wave  characteristics  was  given  by  Melby 
and  Kobayashi  as: 


5(0  =  5(0  +  0.022 


WO* 
( TP)T 


(r0'25-O  M  tn<t<t 


.0.25 


n+ 1 


(3) 


1/2 

where  Nmo  =  Hmo  /  (A DnSo),  Hmo  =  4 (m0)  ,  m0  is  the  zero  moment  of  the  incident  wave  spectrum, 
and  Tp  is  the  spectral  peak  period.  The  empirical  coefficients  in  Equations  2  and  3  will  be 
primarily  a  function  of  structure  slope,  wave  period,  beach  slope,  and  structure  permeability. 
These  equations  have  been  verified  for  the  following  range  of  laboratory  conditions: 


Structure  slope,  tan  a:  1V:2H 

Significant  wave  height,  Hs :  5.05  cm  -  15.80  cm 

Toe  depth,  A:  11.9 -15.8  cm 

Mean  wave  period,  Tm:  1.23  sec  -  1.80  sec 

Iribarren  parameter:  tan  a/(Hs/Lomf'5:  2.08  -  4. 17 

Beach  slope:  1V:20H 

Structure  crest  height,  hc:  30.5  cm 

Stone  density,  pa:  2.66  g/cm 

Armor  stone  gradations:  0.75M50  <  M50  <  1.25 M50,  D85/Di5  =  1.25  and  0.125M50  <  M50  <  4 M50, 
D%5/D\5  =  1.53 

Filter  layer.  (AT5o)armor/(AT5o)filter  —  25,  (f9;)5o)arm0r/(f9;j5o)filter  —  2.9 

The  deepwater  wavelength  computed  from  the  mean  period  is  Lom  =  gT„r/2n,  where  g  is  the 
acceleration  of  gravity. 

Equations  2  and  3  are  plotted  in  Figure  3,  where  S(t )  is  plotted  as  a  function  of  number  of  waves, 
Nw.  Here  Nw  =  t/Tm.  In  Figure  3,  the  measured  mean  eroded  area  plus  and  minus  one  standard 
deviation  are  plotted.  These  data  are  from  one  long  series  of  six  different  storms  with  two  water 
levels.  Equations  2  and  3  should  be  conservative  for  most  applications  because  they  are  based  on 
severely  breaking  waves,  a  relatively  steep  beach  slope,  and  a  relatively  impermeable  core. 
Deviations  from  the  range  of  tested  conditions  are  likely  to  produce  different  empirical  coeffi¬ 
cients  in  Equations  2  and  3.  Caution  should  be  exercised  in  applications  of  these  equations 
outside  of  the  range  of  conditions  tested. 


The  mean  parameters  S,E  ,C,  and  L  and  the  standard  deviations  os,  oe,  oc,  and  <3l  were  used  to 
describe  the  tendencies,  variabilities,  and  ranges  of  damage  and  the  damaged  profile.  All 
measured  values  from  all  measured  series  were  in  the  following  ranges: 

Damage:  -2.7<(S-S)/os<3  (4) 


Eroded  depth:  -2.1<(E-E)/cte<2.1 


(5) 


Cover  depth:  -2.7<(C-C)/crc<2.8 


(6) 
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Figure  3.  Mean  damage  as  a  function  of  number  of  waves  at  mean  period 


These  ranges  allow  the  lower  and  upper  limits  of  the  damaged  profile  descriptors  to  be  estimated. 
In  order  to  reduce  the  number  of  parameters  for  design,  Melby  and  Kobayashi  (1998a)  expressed 
the  key  profile  parameters  as  a  function  of  the  mean  damage  as  follows: 


05=0.55' 0,65 

(7) 

£■=0.465  05 

(8) 

C=Q  — 0.15 

(9) 

Z  =0.445  05 

(10) 

Using  Equations  4  and  7  with  5  =  13,  corresponding  to  localized  failure  in  one  series,  <Js  =  2.65 
and  6  <  S  <  21.  This  illustrates  the  large  alongshore  variability  of  damage  at  failure  for  along¬ 
shore  unifonn  waves. 

EXTENSION  OF  PREVIOUS  MODEL:  A  modification  to  Equation  2  was  introduced  by 
Melby  and  Kobayashi  (1999)  based  on  the  hypothesis  of  equivalency.  The  revised  equation 
allows  for  non-zero  initial  damage  values.  The  modified  equation  is: 
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S(t)  =  0.01  IN s  (N e  +  81V)0'25 


N 


e 


Htn) 

0.01  IN  5 


4 


(ID 


bN  =  (t-  tn  )  /  Tm  for  tn  <  t  <  tn+l 

This  equation  is  similar  to  Equation  2,  but  predicted  damage  is  not  dependent  on  the  time  that  the 
simulation  begins.  Therefore,  future  damage  can  be  predicted  for  an  existing  structure  whose 
history  is  unknown.  This  equation  would  be  appropriate  for  a  major  rehabilitation  study. 

EXAMPLE:  A  single  stonn  example  is  provided  to  show  how  parameters  are  used  in  the 
equations. 

A  traditional  rubble-mound  breakwater  is  to  be  constructed  in  seawater  so  that  the  longitudinal 
axis  of  the  structure  is  parallel  to  the  wave  crests.  The  significant  wave  height  used  for  design  is 
determined  to  be  depth  limited  from  a  shoaling/refraction/diffraction  study.  The  design  wave 
height  is  calculated  at  a  distance  of  5 Hs  seaward  from  the  structure  toe. 

Given  design  values: 

Significant  wave  height:  Hs  =  2.07  m  (6.8  ft) 

Specific  gravity  of  seawater:  S  =  1.0256 
Specific  gravity  of  annor  stone:  S,-  =  2.65 

Specific  weight  of  armor  stone:  yr  =  2.66  x  104  N/m3  (169.6  lb/ft3) 

Density  of  armor  stone:  pr  =  2.72  t/nr  (5.27  slug/ft3) 

Mean  wave  period:  Tm  =  10.8  sec 
Structure  seaward  slope:  tan  a  =  0.5 
Median  stone  size:  Ws o  =  1,31 1  kg  (2,889  lb) 

Stone  gradation:  0.75  W50  <  W50  <  1 .25 IF50 

Storm  duration:  4  hr  =  14,400  sec  (1,333  waves  at  Tm) 

Annor  layer  thickness:  2D„50 

Calculations: 


D 


«50  — 


^50  | 

1/3 

^50 

1/3 

1.311 1 

.  lr  , 

,  Pr  , 

2.72  t/m3 

\ 

0.784  m 


Hs  _  2.07  m 

A Dn50  ~  (2.65-1)0.784  m 


Damage  due  to  a  single  stonn  would  be  computed  using  Equation  2  as 
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S(t)  =  Sft,)  +  0.025^]%  (fa2S -+5) 

(T,„t 

=  0  +  0.25  - : — —  (14,400  s)0'25  =  1.58 

(10.8  s)0'25 

This  damage  level  indicates  that,  for  the  single  4-hr  storm,  there  would  be  an  eroded  area  of 
about  15  to  16  median-sized  stones  removed  from  a  typical  cross  section.  This  is  minimal  damage 
and  would  not  affect  the  integrity  of  the  structure.  Assuming  uniform  wave  height  alongshore, 
the  alongshore  variability,  given  as  the  standard  deviation,  of  damage  would  be  given  by 
Equation  7. 

=  0.5  S  °'65  =  0.5(1 ,58)0-65  =  0.67 

And  the  alongshore  range  of  damage  would  be  given  by  Equation  4  as 
—2. 7(3.0)  +  15.8  <S<  3. 0(3.0)  +  1.58  S>0 

0<S<  3.61 

So  the  range  of  damage  along  the  shore  is  quite  large  but  still  represents  only  minimal  damage  to 
the  structure. 

The  damaged  armor  layer  profile  shape  can  be  analyzed  in  a  similar  fashion.  The  initial  layer 
thickness  is 

tr  =  2 Dn50  =  2(0.784  m)  =  1.57  m 

Then  the  means  of  the  maximum  eroded  depth,  minimum  remaining  cover  depth,  and  maximum 
eroded  cross-shore  length,  respectively,  would  be  given  by  Equations  8,  9,  and  10  as  follows. 

E  =  0.46  S  0  5=  0.46(1. 58)°'5  =  0.58 

C  =  Ca  -  0.1  S  =  1.57  -  0.1(1.58)  =  1.41 

L  =  0.44  S  0  5  =  4.4(1.58)°  5  =  5.5 

The  mean  dimensional  eroded  depth  is  de  =  E  Dn 5o  =  0.58  (0.784  m)  =  0.45  m.  The  mean 
dimensional  minimum  remaining  cover  depth  is  dc  =  CD, ,5 o  =  1.40  (0.784  m)  =  1.11  m,  so  there 
is  significant  cover  over  the  underlayer.  Note  that  the  initial  thickness  of  the  armor  layer  is  not 
uniform;  so  the  maximum  eroded  depth  may  not  occur  in  the  same  location  as  the  location  of 
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minimum  remaining  cover.  In  addition,  the  sum  of  E  and  C  may  not  be  equal  to  Ca  for  the  same 

reason.  The  mean  of  the  maximum  eroded  length  or  cross-shore  eroded  hole  is  le  =  4.34  m.  This 
stone  size  was  very  well  suited  for  this  application  if  this  was  a  design  level  stonn. 

If  a  storm  of  Nw  =  5,000  waves  were  to  strike  the  undamaged  structure,  then  t  =  5,000  (10.8  s) 
=  1 5  hr.  Then  the  mean  damage  would  increase  to 

_  i  za5 

S(t )  =  0  +  0.025  -  (54,400  s)0'25  =  2.20 

(10.8  s)0'25 

This  damage  level  is  still  quite  acceptable  for  most  situations. 

CONCLUSIONS:  This  CHETN  provides  a  method  for  computing  damage  on  a  rubble-mound 
breakwater,  revetment,  or  jetty  trunk  section  stone-armor  layer  that  is  exposed  to  depth-limited 
breaking  waves.  The  methods  are  useful  in  determining  the  expected  life  of  new  or  damaged 
structures  and  in  developing  life-cycle  analyses.  The  empirical  equations  should  be  conservative 
for  most  applications,  but  care  should  be  exercised  when  applying  outside  of  the  tested 
conditions,  as  described  in  this  CHETN. 

POINTS  OF  CONTACT:  The  author  would  like  to  thank  Dr.  Shamsul  Chowdhury  for 
valuable  technical  assistance  with  editing  this  document.  This  CHETN  was  produced  with 
funding  from  the  Navigation  Systems  R&D  Program  within  the  Prediction  and  Prevention  of 
Breakwater  Deterioration  Work  Unit.  The  Principle  Investigator  for  the  Work  Unit  was 
Dr.  Jeffrey  A.  Melby.  Questions  about  this  CHETN  can  be  addressed  to  him  at  (601-634-2062  or 
email  Jeffrev.A.Melbv@erdc.usace.armv.mil). 
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